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ABSTRACT 

Excess emission from a point-like source coincident with the central star of the Helix Nebula is 
detected with Spitzer at 8, 24, and 70 /im. At 24 ^m, the central source is superposed on an extended 
diffuse emission region. While the [O IV] 25.89 /Ltm line contributes to the diffuse emission, a 10- 
35 /zm spectrum of the central source shows a strong thermal continuum. The excess emission from 
the star most likely originates from a dust disk with blackbody temperatures of 90-130 K. Assuming 
a simple optically thin debris disk model, the dust is distributed in a ring between ^35 and ~150 
AU from the central star, possibly arising from collisions of Kuiper-Belt-like Objects or the break-up 
of comets from an Oort-like cloud that have survived from the post-main-sequence evolution of the 
central star. 

Subject headings: infrared: stars - stars: individual (WD 2226—210) ~ white dwarfs 
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1. INTRODUCTION 

The Helix Nebula (NGC 7293), at only 219 pc 
([Harris et al.ll2007[ ). is one of the best studied planetary 
nebulae. Hubble Space Telescope (HST) images at optical 
wavelengths and recent Spitzer Space Telescope (Spitzer) 
images at 3.6-8 /xm have resolved the Helix Nebula into 
thousands of "cometary knots" immersed in diffuse gas, 
produ cing strikingly be autiful and w ell-publicized images 
(jO'DeU et al. 2004: H ora et al.ll200^ . In contrast, Chan- 
dra X-ray observations of the Helix detected only a point 
source associated with the central star, but the X-ray 
properties of this source are mysterious. 

The central star of the Helix Nebu la, WD 2226-210 
is a ~110,00q K ho t white dwarf l|Napiwotzkil Il999t 
iTraulsen et al.l l2005l ). As expected, its photosphere 
emits soft X-rays at energies up to ~0.25 keV; how- 
ever, Chandra and ROSAT observations revealed ad- 
ditional hard X-ray emissi on near 1 keV (jLeahv et al.l 
119941 : iGuerrero et ani2001h . The luminosity, temporal 
variations, and spectral propert ies of this hard X-ray 
component (|Guerrero et al.|[2001h and th e variability of 
WD 2226-210's Ha emission line profile (jCruendl et all 
l2001h are consistent with the coronal activity of an M 
dwarf. However, a sensitive HST search has ruled out 
the presence of any companion with a spectral type ear- 
her than M5 (|Ciardullo et al.lll999f ). 

The mystery regarding the central star of the Helix 
Nebula deepens with new Spitzer observations at 24 and 
70 /im, as reported in this paper. 
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The Helix Nebula was observed with all three 
Spitzer inst r umen ts: the Infrared Array Camera (IRAC; 
iFazio et al.l I2004D the Infrared Spectrograph (IRS; 
Hou ck et all I2004D . and the Multiband I maging Pho- 
tometer for Spitzer fMIPS: [Rieke et al.ll2004f) . The MIPS 
observations were made in the scan map mode with a 
medium scan rate and a half-array offset, resulting in 
an effective integration of 80 s per pointing and cov- 
ering an area of 0?5 x 0?8 in each of the 24, 70, and 
160 fim bands. Three MIPS scan maps were taken in 
sequence with >10 hr separations to avoid confusion by 
asteroids. The IRS high-resolution observations of the 
central star and a nearby nebular position were made in 
sequence in the staring mode. The IRAC observ ations 
have been described in detail bv iHora 

Basic reduction for the MIPS dat a was carried o ut us- 
ing the MIPS Data Analysis Tool (jCordon et al.lf200a . 
IRAC mosaics were generated from the Spitzer Science 
Center (SSC) Basic Calibrated Data (BCD) product 
(ver. SI 3. 2) using a customized IDL program. The 
central 15.'4 xl2'.5 region of the nebula at 3.6, 4.5, 
8.0, 24, 70 and 160 /xm is displayed in Figure [T] IRS 
high-resolution spectra were extracted from the SSC 
BCD product (ver. S14.0) using the IRS team's Spec- 
tral_ Modeling, Analysis, and Reduction Tool (SMART; 
Higdon et al. 2004). As the spectra on the star were 
extracted using a full-slit width mode in both the Short- 
High (SH, 9.9-19.6 /Ltm, 5"x 11" slit) and Long-High (LH, 
18.7-37.2 iim, ll"x22" slit) modules without off-source 
sky subtraction, the LH spectrum was scaled to match 
the continuum level of the SH spectrum. The scaled LH 
spectrum and the SH spectrum are plotted together in 
the top panel of Figure [2] 

The IRAC images of the Helix Nebula show a bare 
point source at the location of the central star, WD 
2226—210 (see Fig. 1). This central source is de- 
tected with flux densities of 374±19, 241±24, 171±26, 
and 174±17 /xJy at the 3.6, 4.5, 5.8, and 8.0 /im 
IRAC bands, respectively. The flux densities in the 
first three bands are, within error limits, consistent 
with those expected from the Raylcigh- Jeans tail of a 
110,000 K blackbody model normalized to the Two Mi- 
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Fig. 2. — Spitzer spectrum of the Helix central star. The top 
panel shows the spectrum without sky (nebula) subtraction; there- 
fore, the spectrum contains contributions from nebula, zodiacal 
background and the central point source. For comparison, the 
MIPS 24 /im bandpass is shown as a dotted line in the top panel. 
The bottom panel shows the spectrum from the central souce af- 
ter nebular and sky subtraction and Gaussian smoothing (with no 
scaling between SH and LH modules). MIPS 24 /im photometry is 
shown as a square. A blackbody of 102 K is superposed as a long 
dashed line. 
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Fig. 3. — SED of the central white dwarf in the Helix Nebula. 
The photometry for ground-based Johnson B and V l IHarris et al.l 
l2007l ) is displayed as open circles, as filled diamonds for 2MASS, 
filled circles for IRAC, and open diamonds for MIPS photometry 
with a downward arrow showing the 160 fj,m 3 cr upper limit. The 
white dwarf photospheric model is shown as a grey line. The ex- 
cess SED is consistent with a blackbody with temperatures (BBT) 
ranging from 86 K (the best fit using MIPS 24 and 70 ^m points) 
to 128 K (the best fit using IRAC 8 /^m and MIPS 24 fira). An 
intermediate temperature of 102 K is also shown as a heavy dotted 
line for comparison. The SED of a dusty ring (35-150 AU) is shown 
as a (red) long-dashed line. 



cron All Sky Survey (2MASS) photometry of the star 
(2MASS 22293854-2050136), J = 14.332±0.024, H = 
14.490±0.039, and = 14.548±0.079. The flux density 
at 8 /im is, however, twice as high as expected. 

The MIPS images of the Helix Nebula show very differ- 
ent distributions of emission. In the 24 band, a point- 
like source coincident with the central star is superposed 
on a plateau of diffuse emission. In the 70 /im band, 
no diffuse emission is seen, and only a central pointlike 
source is detected. In the 160 /im band, neither the cen- 
tral source nor the bright diffuse emission is detected. 

To analyze the central point-like source at 24 /^m, the 
bright diffuse background needs to be subtracted. The 
surface brightness of the diffuse emission falls off linearly 
with radius (r) from the central star. Therefore, we use 
a linear model fit to the surface brightness profile be- 
tween r — 31" and 100" to estimate the diffuse emission 
superposed on the central point-like source. After back- 
ground subtraction, the central source appears slightly 
resolved with a FWHM of ^--^9", 1.5 times that of a true 
point source, possibly due to an imperfect subtraction 
of the background. Using aperture photometry and a 
color correction for a ~100 K blackbody (see H3.3\ for 
justification), the background-subtracted flux density of 
the central source is 48.4±7.3 mJy in the 24 fim band. 
No variability is detected at 24 /xm in our three observa- 
tions. Using a similar method, we measure a flux density 
of 224±33 mJy in the 70 /im band and a 3 ct upper limit 
of 237 mJy in the 160 fj,m band. 

3. DISCUSSION 

3.1. The Diffuse Emission and Point-Source Regions 

The detection of a point-like source in the MIPS 24 and 
70 /im bands is puzzling because the large extent of the 
ionized nebula is well resolved with Spitzer. While the 



MIPS 24 /im flux from the unresolved source could origi- 
nate from emission lines in the bandpass, the lines would 
have to originate in high-ionization stages from abundant 
elements. Possible lines are the [Ne V] 24.32 /im line and 
[O IV] 25.89 /xm line , the l atter having been suggested 
W iLeene Sz PottaschI (|1987f ) from the /i?yl5 observation. 
The ionization thresholds needed to produce O IV (54.9 
eV) and to excite the A4686 line by the ionization and 
recombination of He II (54.4 eV) are similar. Therefore, 
we resort to a He II A4686 image of the Helix Nebula 
(jO'DeU et all [200I to assess the spatial extent of the 
[O IV] emission at a higher resolution. The He II image 
of the Helix Nebula shows a central diffuse emission re- 
gion similar to that seen at 24 /im, but not a point-like 
central source. From this comparison we conclude that 
the [O IV] line contributes to the diffuse emission up to 
a radius of ~2'. 

The absence of the central point-like source in the He II 
image confirms that the [O IV] line does not contribute 
to the central peak in the MIPS 24 /im band. Another 
promising candidate is the [Ne V] 24.32 /im line. With 
an excitation potential of 97 eV, Ne V can be produced 
via photoionization by hard X-rays or thermal coUisional 
ionization at temperatures (T) higher than 10^ K. At 
such high temperatures, the He II A4686 recombination 
line may be weakened to below the detection limit be- 
cause of its r~i/^ dependence. While this might explain 
the absence of the central peak in the He II image, it 
cannot explain the central point-like source in the MIPS 
70 /im band, as this band does not contain any nebular 
lines of abundant elements at high ionization stages. 

Both [O IV] and [Ne V] lines are seen in our IRS spec- 
trum of the central region (see the top panel of Fig. [2]) 
on top of a continuum (a combination of dust grains in 
the ionized nebula, zodiacal background and the central 
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point source). Without a proper subtraction of the neb- 
ular and zodiacal contribution, it is difficult to assess the 
relative importance of the emission lines and the contin- 
uum from the central star based on this IRS spectrum 
alone. The true spectrum of the central source is re- 
vealed after subtraction of the nebular and sky contribu- 
tion (shown in the bottom panel of Fig. [J). The spectrum 
of the central source will be addressed in more detail in 
another paper. 

3.2. Nature of the Central Infrared Source 

The spectral energy distribution (SED) of the cen- 
tral source is shown in Figure [3l For the white dwarf, 
we employ TMA P, the Tubingen NLTE Model Atmo- 
spher e Package (jWerner et al.l 120031 iRauch fc Deetjenl 
|2003D , in its most recent version for the calculation of 
plane-parallel, static white dwarf models that consider 
H-He-C-N-O-Mg-Si-Ca-Ni model atoms. The SED 
shortward of 6 /xm is consistent with the scaled white 
dwarf photospheric model for 110,000 K. Assuming that 
the excess infrared emission originates from a blackbody 
radiator, its temperature is 90-130 K (see Fig. [3|), far 
too low for a stellar object. If the blackbody radiator is 
heated by the central star WD 2226—210, of which the 
effective temperature is 110,000 K and radius is 0.024 
i?0, it would have to be 41-91 AU from the star. This 
radial distance corresponds to 0'.'19-0'.'42, so it cannot be 
resolved with Spitzer. The luminosity of this blackbody 
radiator, (5.1-11) xlO^^ ergs s~^, requires a total cross 
section (atot) of (0.8-8.7) xlO^^ cm^, or 3.8-38 AU^. 
Although we do not have any constraint on the geome- 
try of the emission source, only an extended dust cloud, 
most probably a disk can provide such a large emitting 
area. The extinction toward the star is smal l based on 
the International Ultraviolet Explorer data (Bohli n et al.l 
[1982) ■ consistent with the 219 pc distance. The dust disk 
is likely to be optically thin and tilted since the main ncb- 
ular ring is til ted by ~35° with respect to the line of sight 
(| Young et al.l fr999). 

Given the evolution of the Helix central star through 
the asymptotic giant branch (AGE) and planetary neb- 
ula phases, the presence of a debris disk from the 
stellar main-sequence phase is surprising but not en- 
tirely impossib le. Recently, dust disks around the w hite 
dwarf GD 362 (jBecklin et al .1'2005': 'Kilic et al."2005') and 
the white dwarf G29-38 fZuckerman & Bccklin 1987; 
iReach et al.ll2005l ) were reported; these disks are closer 
to their central stars and thus have SEDs showing excess 
emission starting from the near-IR bands. The Helix cen- 
tral star is much hotter and more luminous than these 
two white dwarfs (i.e., the Helix central star is much less 
evolved along its cooling track than the other two) ; thus 
it can heat a debris disk at a larger radial distance. 

3.3. Physical Properties of the Debris Disk 

If the infrared excess around the Helix central star in- 
deed arises from a debris disk, the dust grains in the disk 
that dominate the radiation are likely to be microns to 
millimeters in size. The minimum grain size will be set 
by blowout by radiation pressure . Given t he stellar mass 
of 0.58 Mq (Napiwotzki 1999; Tr aulsen et al. 2005) and a 
luminosity of 76 Lq (from the observed photometry and 
distance), grains smaller than ~60 fim will be ejected 
from the system. For an initial estimate, one can assume 



that the disk consists of grains with sizes (radius of a) 
from 60 to 1000 fim and has a power-law size distribu- 
tion, n(a) PC Q~^'^i a s expec ted for a coU ision-dominated 
disk ()Dolinanvilll"969 : Tanak a et a"nil996f ). The relatively 
large grains support using blackbody-based approxima- 
tions to characterize the system. 

For a simple estimate, we take the disk dimensions 
from the preceding section and approximates the grain 
size distribution by a single size of 100 ^m, character- 
istic of the grains that dominate the infrared emission. 
The dust mass in these grains is (4/3) PgOatot where 
Pg is the grain density and atot is the total area that 
intercepts the star light. Since the total dust mass is 
proportional to J n{a) a^da oc a"-^, the total dust mass 
from the disk is '^O.ll Af^ assuming pg=2.5 g cm~^ and 
ftot =6xl0^''cm^. 

This rough estimate can be tested with a more real- 
istic model that fits the SED, although such models are 
still subject to degeneracy and ambiguiti es. We assume 
that the grains are astronomical silicates (jLaor fc Draind 
[1991 with a size distribution as before. The disk is op- 
tically thin, is heated only by WD 2226-210, and has a 
constant surface number density, with inner and outer 
boundaries of ri„ and rout- The observed SED is fitted 
by a ring extending from 35 to 150 AU with a total dust 
mass of O.13Af0(±5O%), as shown in Figure [3l 

The exact inner radius of the disk is subject to the 
grain properties used in the model; however, the inner 
radius cannot be too close to the star because no observ- 
able excess is found at wavelengths shorter than 8 pm. 
An inner hole is required by the SED modeling. The 
outer boundary of the disk is constrained in our model 
because a disk with a much larger outer radius would 
have higher 160 pm flux, contradicting the observed up- 
per limit. Therefore, the dust around the Helix central 
star is in a form of ring between ~35 and ~150 AU. This 
configuration is remarkably similar to the planetary de- 
bris disks seen around main-sequence stars. 

3.4. Origin of the Debris Disk 

What is the origin of the dust Spitzer detects around 
the Helix central star? It has been established that any 
planet closer t han ^1 AU will be en gulfed by an expand- 
ing red giant (Sies s fc Liviolll99"9al [ai. while planets out- 
side '^5 AU from the Sun wil l surviv e post-ma in-sequence 
evolu tion (jSackmann et al.l Il993t [Pebes fc SigurdssonI 
I2002D . and the orbits of surviving planets and most of 
the Kuiper Belt objects (KBOs) and Oort Cloud comets 
will expand adiabaticall y and remain bound to the so- 
lar system (iDuncan fc L issaucr 1998). The re-stabilized 
KBOs and Oort Cloud can later become the source of 
objects that go into the inner part of the system, either 
plunging into the white dwarf or breaking up due to tidal 
destruction, and they can populate the inner system with 
dust. A similar scenario has been proposed to explain the 
photospheric meta l contamination in DAZ white dwarfs 
(IStern et all [l990l : IParriott fc Alcoca [19981: [Jural [2006h 
and may even prod uce the observed X-ray eni ission at 
1 keV. Furthermore. iDebes fc SigurdssonI (|2002D suggest 
that if the surviving planet system becomes unstable af- 
ter the AGB phase, the entire system becomes dynami- 
cally young; i.e., frequent collisions and encounters can 
occur among surviving planets and comets, leading to a 
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period of enhanced "late bombardment" . It is likely that 
the kilometer-size objects surviving around the Helix cen- 
tral star are experiencing similar perturbations and pro- 
ducing fine dust in the system. The Poynting-Robertson 
(P-R) hfctime is 3 - 50 x 10^ yr for -100 ^im grains lo- 
cated at 35 and 150 AU, while the coUisional lifetime is 
0.1 — 1 X 10^ yr. Note that these timescales are much 
lon ger than the kin ematic age of the nebula (2.8x10^ 
yr; lNapiwotzkilll999f ) and a typical timescale in the post- 
AGB evolution (— 4 x 10* yr for a 0.565 Mq white dwarf; 
lSchoenberneijll983[ ). Assuming an average dust lifetime 
of 1 X 10^ yr, the dust production rate is 2.5 x 10^'^ g s~^. 
This rate is only a lower limit since the dust production 
is not in equilibrium, but it is ~10^ times the dust pro- 
duction rate observed from comet Hale-Bopp, consistent 
with the hypothesis of an enhanced late bombardment 
in the post-main-sequence evolution. 

The total observed dust mass, ^^0.13 Mq, is equiva- 
lent to 20 asteroidal objects with a radius of ~1200 km. 
While the dust cloud around the white dwarf G29-38 
may arise from tidal disruption of Pluto-sized objects 
((jura 2003), this cannot be the case for the disk in WD 
2226-210 because the tidal disru ption zone is much closer 
to the central star, -0.003 AU i| Jural l2003l . eq.[4]). It is 
thus unlikely that the debris disk is produced by tidal 
breaking of large asteroidal objects, but rather by en- 
hanced collisions among the surviving Kuiper Belt-like 
and/or cometary objects. 

The inner radius of a P-R drag-dominated disk, in gen- 
eral, should be very close to the star, up to the dust sub- 
limation radius (~0.1 AU in the case of the Helix central 
star) if there is no planet-size or sub-stellar object inside 
the disk to stop the dust grains from spiraling inward. 
The inner hole (r<35 AU) derived from our SED model 
is likely caused by the dynamical perturbation of un- 
seen planets or sub-stellar objects inside the hole, which 
can trap the migra ting dust or scatter the dust outward 
(|Mever et al.ll20d^ . Alternatively, the inner hole may 
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also be explained by ice sublimation if the grains are icy 
(jJura et al.l Il998^ ■ a process that will occur at a dust 
temperature of ~150 K, corresponding to ~30 AU. 

If the Helix central star had a close binary in the main- 
sequence phase, a binary merger event could produce an 
accretion disk around the remaining white dwarf. The 
physical structure of the disk would be highly dependent 
on the properties of the original binary system; however, 
the fact that there is no near-infrared excess suggests 
that the inner part of the optically thick accretion disk 
would have to dissipate in a short period of time after the 
binary merged. It is hard to explain why the dust is so 
far away from the white dwarf if it was the remnant disk 
due to close binary evolution. Alternatively if the com- 
panion was in an intermediate orbit and avoided spiraling 
into the primary during the evolution, the binary config- 
uration is efficient at trapping the dust from the mass 
loss of the primary t o form a stable circumbinary disk in 
the p ost-AGB phase (jvan Wincke]|l2003l : Ide Ruvter et aH 
120061 ). inducing the presence of the dusty disk we de- 
tected. However, any of the micron-size dust generated 
in the AGB phase should be blown out by the strong 
radiation in the early white dwarf evolution. The fact 
that we do see dust around the Helix central star favors 
the debris disk scenario. Future high spatial resolution 
infrared imaging is needed to further constrain the ge- 
ometry of the emitting region. 
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Fig. 1. — Spitzer images of the Helix Nebula. Each panel is displayed in a false-color logarithmic scale with brightness and contrast 
adjusted for best presentation and oriented as north up and east toward the left with a field of view of 15f4xl2.'5. 



